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The path integral formahsm is apphed to derive the full partition function of a generalized Su- 
Schrieffer-Heeger Hamiltonian describing a particle motion in a bath of oscillators. The electronic 
correlations are computed versus temperature for some choices of oscillators energies. We study the 
perturbing effect of a time averaged particle path on the phonon subsystem deriving the relevant 
temperature dependent cumulant corrections to the harmonic partition function and free energy. 
The method has been applied to compute the total heat capacity up to room temeperature: a low 
temperature upturn in the heat capacity over temperature ratio points to a glassy like behavior 
ascribable to a time dependent electronic hopping with variable range in the linear chain. 
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1. INTRODUCTION 

Determining the polaron properties is a relevant problem in many body theory [1,2]. The electron motion through 
the crystal is generally accompanied by a lattice deformation whose size and entity depend both on the strength of 
the electron-phonon coupling and on the value of the adiabaticity parameter peculiar of the system [3-18]. In general, 
the lattice deformation does not follow instantaneously the electron dragging it and the retardation effect becomes a 
key ingredient of this many body problem. The Holstein [19] and the Su-Schrieffer-Heeger (SSH) [20] Hamiltonian 
models are fundamental tools in polaron physics. While the former has been mainly used to describe electron coupling 
to local optical phonons, both acoustic [21] and optical branches are present in the dimerized lattice of the latter. 
However, acoustic phonons in the SSH model arc not essentially affected by the e-ph coupling whereas optical phonons 
are strongly softened due to single electron polarization [22] . 

In this paper we address the general problem of the interplay between an electron and the optical phonon subsystem 
(a set of independent oscillators providing the bath for the particle) in a generalized SSH model, studying i) the 
temperature dependent electronic correlations induced by the dissipative phonon bath and ii) the thermodynaniical 
behavior of the latter in the presence of a perturbating particle motion. The path integral method [23] , being valid 
for any value of the e-ph coupling, seems particularly suitable to our task. Moreover, it naturally introduces the time 
(as an inverse temperature) into the formalism thus allowing us to derive the thermal particle correlation functions 
and the cumulant expansion for the phonon free energy. In Section 2, we present the model and discuss the particle 
correlations in the oscillators bath. The phonon free energy is calculated in Section 3 while in Section 4 we apply 
the model to compute the total heat capacity of the system both in weak and strong e-ph coupling regimes. Some 
conclusions are drawn in Section 5. 

2. THE MODEL 

The ID SSH interacting Hamiltonian reads: 

r 

Jr,r+1 = + ^ "r+l)] (0-1) 

where J is the hopping integral for an undistorted chain, a is the electron-phonon coupling, is the dimerization 
coordinate relative to the displacement of the atomic group on the r— lattice site along the molecular axis, /J and 
fr create and destroy electrons on the r— group. The free Hamiltonian is given by a set of classical independent 
oscillators. By introducing x{t) and ?/(t') as the electron coordinates at the r and r + 1 lattice sites, respectively and 
mapping Ur — > u{t) and ii(''"') we transform the real space Hamiltonian of eq.(l) into the time dependent 

Hamiltonian: 
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H{t,t') = Jr,r'[f^{x{T))f{y{T')) + f\y{T'))f{x{T))) 



-^ry = --[J + aiuiT)-u{T'))] 



(0.2) 



The SSH Hamiltonian has a twofold degenerate ground state which undergoes a Peierls instabihty. In real space 
the soliton connects the two degenerate phases with different senses of dimerizations and a localized electronic state is 
associated with each soliton. Both electron hopping to band states (thermal excitation) and electron hopping between 
solitons are allowed. 

Mapping the Hamiltonian onto the time scale we set up the finite T formalism in which thermally activated electron 
hops become time dependent and the retarded nature of the interactions is accounted for. H{t, t') is more general 
than the real space SSH Hamiltonian since hopping processes are not constrained to first neighbors sites along the 
chain. A variable range hopping introduces some local disorder in the system [24] , therefore the present approach may 
apply to polymers in which hopping conduction mechanisms prevail. We point out that in real systems, the Peierls 
gap can be smeared by temperature or doping induced disordering effects [25]. Eq.(2) displays the semiclassical nature 
of the model as quantum mechanical degrees of freedom interact with the classical variables u{t). Setting t' — 0, 
w(0) = 2/(0) = 0, averaging the electron operators over the ground state and pinning the chemical potential to the 
zero energy level, we write the average energy per lattice site: 



< H{t) > 



N 

V{x{t)) = -J 



= V{x{t)) +w(r)i(r) 

n/a 



dk cos[fca;(T)] cosh.{ekT/h)nF{ek) 



j{T) = —a— I dkcos[kx{T)]cosh{ekT/fi)nF{ek) 
7^ Jo 



(0.3) 



where N = L/a, with L the chain length and a the lattice constant, np is the Fermi function and Ck = — Jcos(fc) 

is the electron dispersion relation, ^(.^(t)) is an effective term accounting for the t dependent electronic hopping 
while jir) is the external source current for the oscillator path u{t). Averaging the electrons over the ground state we 
neglect the fermion-fermion correlations [26] which lead to effective polaron-polaron interactions in non perturbative 
analysis of the model. This approximation however is not expected to affect substantially our thermodynamical 
calculations. Being the energy in eq.(2) linear in the displacements we can write the general path integral [27] at any 
temperature as: 



< xiPMo) >= n/ ^"i(^) / ^^(^) ■ 

\fj dT{^xHr) + V{x{T))-Y,u,{T)3{T] 



exp 



exp 
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where wc have taken a large number of oscillators (^^(t), i = I..N) as the bath for the quantum mechanical particle 
whose coordinate is x(t). (3 is the inverse temperature, m is the electron mass and are the oscillators frequencies. 
The oscillator masses are considered as independent oii, Mi = M and hereafter we set M = lO^m. After integrating 
out the oscillators coordinates over the paths Duiir), imposing a closure condition ^a;(/3) = a;(0)^ on the particle 
paths and replacing t ^ t/U, we obtain the full partition function in the functional form 



■^(iW) = Zph j) Dx{T)exp 



.-x\r)-V{x{r))-A{j{r)) 



N 



Zph = JJ 
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The thermodynamics of the interacting system can be derived from eqs.(5) as a function of J, a and oscillators 
bath. An application is shown in Section 4 where the total heat capacity is computed versus temperature. In terms 
of the generating functional Z{j{T)) the two particle correlation function is defined as 



(0.6) 

Ji=o 



Then, using cqs.(5) and the fact that the source action A{j{T)) is quadratic in the current, we can study the effect 
of the phonon bath on the electronic time correlations at any temperature. Setting T2 = 0, we plot in Figures 1 the 
square root of G^^\t) =< x{t)x{0) > (r e [0,/3]) for three choices of oscillator baths: i) a low phonon spectrum 
made of ten oscillators Uuii = 2meV.....fiu)io = 'iQmeV (spaced by 2meV), ii) an intermediate phonon spectrum with 
TiLJi = 22meV.....fiLUia = AQmeV and iii) an high phonon spectrum with fiuji = A2meV.....huJif) = QOmeV . N and 
oscillator energies are purely arbitrary. The periodic condition G^^^ (t) = G^^^ (/3 — t) holds as a general property of the 
particle-particle correlation function which, up to a prefactor, coincides with the sum of the oscillators Green functions 
and therefore it docs not depend on the width of the electron band. At low temperatures (Fig. 1(a)), small oscillator 
energies favor larger electronic correlations which however get the maximum value of 0.26A at r ^ 0. Increasing the 
temperature at T = lOOii' (Fig. 1(b)) leads to a substantial enhancement of i/G(^)(t) which becomes of order lA 
in the case of a low energy phonons bath. At room temperature (Fig. 1(c)) and for the same bath the electrons are 
correlated over a distance of about 2A in the whole r-range. Larger phonon spectra tend to lock the electrons thus 
reducing the particle mobility over the time scale and decreasing the space electronic correlations. Since the upper 
limit in the t range of Figs. 1(a)- 1(c) is necessarily T dependent we plot in Fig. 2, for an interm.ediate phonon bath, 
\/G(^)(t) versus a common r-axis to emphasize the temperature effect on the particle correlations. Only the room 
temperature correlation function retains periodicity over the selected r-range. 

While we focus in the following on the equilibrium thermodynamics of the particle-phonons interacting system, 
it should be noted that generalization of eq.(4) through the closed-time path formalism [28] would permit to derive 
dissipative properties due to the phonon bath friction from the two points correlation function. A density matrix 
study for the similar problem of a particle in a photon bath (in three dimensions) has been carried out by Haba and 
Kleinert [29]. 



3. PHONON FREE ENERGY 



Next we turn to analyse the effect of the particle motion on the phonon subsystem. The set of N independent 
oscillators is perturbed by the particle path x{t) which couples to each oscillator through the strength a (assumed 
independent of i) of the SSH Hamiltonian. Then, the partition function of the disturbed oscillators system can be 

expanded in perturbation series as 



Zph[x{T)] ~ (1- < G > + < G^ > - < G^ > + • •)Zi 



< C' >-- 



-TT f or r f'^ 

^ph n f ■^'"'^'^^ fcT J dTiUi{Ti)x{Ti) ■ ■ J dTkUi{Tk)x{Tk) ■ exp 



dT 



r)+a;,V(T)) 
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In general, the total phonon partition function in presence of an ensamble of particle paths [a;(T)] would be given 
by: Z^f^ = f Dx{T)Zph[x{T)]. Here we take a single particle path, approximating it by the averaged (dimensionless) 

value < x{t) >= ^ dTx{T) = xq. Then, the odd k-terms in the series expansion vanish. Since the oscillators 
are decoupled (and anharmonic effects mediated by the particle path are here neglected) we can study the behavior 
of the cumulant terms < G'' > by selecting a single oscillator having energy lo and displacement u{t). Hence, after 
expanding the oscillator path in Np Fourier components 



Nf 

u{t) = Uo + ^2 ^ (^^n cos(a;„T) — sin(ti;„T) 



(0., 



3 



with ujn = Inn/ fi and taking the measure of integration 

--^\ aWf, /'27r • •27VF7r)^ ^oo Nf 



(0.9) 



V^Xj^ "-00 

with Am = we obtain for the k — th cumulant the following expression 

^jv^ P'' A;! 7rfc/2(a;/3)fe+i (27r)2 + (w/3)2 ' ' (2JVp7r)2 + (a;/?)^ 

Since the cumulants should not depend on the number of Fourier components in the particle path expansion, 
eq.(lO) provides a criterion to set the minimum Np (at any temperature and for any oscillator) through the condition 
2Nftt ^ Lo(3. The thermodynamic properties of the perturbed oscillator can be computed observing that the cumulant 
corrections to the free energy are given by: 



F = -— In 

^ cum p '■^^ 



1+ < > + < > + < > +.. 



(0.11) 



In Fig. 3(a), the three lowest order cumulants are plotted for a weak e-ph coupling and an high frequency oscillator. 

At very low temperatures, the cumulants attain the largest values which are substantially independent of the fc-order. 
By increasing T all cumulants decrease, < > yields the dominant contribution and the perturbation expansion 
converges rapidly. As shown in Fig. 3(b), the free energy corrections are not relevant due to the weak a coupling. Note 

that Fcurn includes Fph plus the perturbation given by eq.(ll) to the fc-order. The k = A correction is not appreciable 

with respect to Fcum- 

In Figures 4 we take a sizeable e-ph coupling (whose order of magnitude applies to a polymer as polyacetylene), 
a factor ten larger than in Figs. 3: at low T, high order cumulants have to be included in the series expansion but 
the temperature range (in which their contribution is relevant) shrinks by increasing k. An high number of Fourier 
components is required in order to get numerical convergence (i.e., Np ^ 18000 at T = IK) signalling that the 
particle dynamics strongly interferes with the oscillator at low T. For any a and any T, one determines the fc-order 
which makes the cumulant series convergent. For instance, being at T ~ 150K: < >^< >~< >, the series 
can be truncated for A: = 6 at T > 150K. Instead, at T ~ 90K, < >« »-< > so that the fc = 8 
term suffices in the series expansion at larger temperatures. Note however that high order cumulant terms should be 
handled with care in computation of physical quantities since they may provide some significant contributions to the 
free energy, mainly at growing temperatures. As it is shown in Fig. 4(b), the largest correction to the free energy is 
due to the k = 2 term which strongly reduces the harmonic value in the whole temperature range whereas the A; = 4 
and k = 6 contributions are also relevant although of decreasing importance. The k = 8 correction, still appreciable 
in the range 5QK :< T ^ 20QK, tends to vanish above room temperature. By enhancing the temperature, cumulants 
and free energy become numerically stable taking a smaller Np'- at T = 300iir, Np ~ 60. 



4. HEAT CAPACITY 



By mapping the electronic hopping motion onto the time scale, we have introduced a continuum version of the 
interacting SSH Hamiltonian . Unlike previous [24] approaches however, our path integral method is not constrained 
to the weak e-ph coupling regime and it can be applied to any range of physical parameters. To compute eqs.(5) 
one has to select the class of particle paths which mainly contribute to the partition function and fix the physical 
quantities characterizing the system: the bare hopping integral J, the oscillator frequencies u>i and the effective 
coupling X = ci^h^/M (in units meV^). 

We take here a narrow band system (J = lOOmeF) to be consistent with previous investigations [30] and with the 
caveat that electron-electron correlations may become relevant in narrow bands. The total heat capacity has been first 
(Fig. 5) computed up to room temperature assuming the low phonon spectrum of Fig.l. The lowest energy oscillator 
yields the largest contribution to the phonon partition function mainly in the low temperature regime while the loiq 
oscillator essentially sets the phonon energy scale which determines the size of the e-ph coupling. A larger number N 
of oscillators in the aforegiven range would not significantly modify the calculation. 

In the discrete SSH model, the value a = /{ttkJ) ~ 1, marks the crossover between weak and strong e-ph 
coupling, with k being the effective spring constant. In the continuum and semiclassical model of eq.(5) the effective 
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coupling is the above defined x- Although in principle, discrete and continuum models may feature non coincident 
crossover parameters, wc assume that the relation between a and J obtained by the discrete model crossover condition 
still holds in our model. Hence, at the crossover we get: Xc ^ ttJU^uj^q/QA. This means that, in Fig. 5, the crossover is 
set at Xc ~ 20QQmeV^ . The total heat capacity over temperature ratio shows a peculiar low temperature upturn (also 
in the weak x regime) which can be mainly ascribed to the sizable effective hopping integral term V{x{t)) . The e-ph 
coupling however determines the shape of the low T anomaly. The small phonon contribution to the heat capacity is 
also reported on to point out that the Dulong Petit value is achieved at T ^ 2007^". 

The effect of the oscillators on the heat capacity is pointed out in Fig. 6 where we take the intermediate phonon 
bath with energies: fioji = 22m,eV, ...,huJio = AOm,eV. Accordingly the crossover is set at Xc ^ SOOOmeV^ and three 
plots out of five lie in the strong e-ph coupling regime. The heat capacity grows fast versus temperature at strong 
couplings due to the source action contribution whereas the presence of the low T upturn in the total heat capacity 
over T ratio is confirmed. Note that, due to the enhanced oscillators energies, the phonon heat capacity saturates at 
T ~ 4007^. 

Integrating eq.(5), we select, at any temperature, the ensamble of particle paths over which the hopping potential 
V[x(t)) is evaluated. This ensamble is therefore T dependent. However, given a single set of path parameters one can 
monitor the V(^x{t)^ behavior versus T. It turns out that the hopping decreases by lowering T but its value remains 
appreciable also at low temperatures 20K). Since the dr integration range is larger at lower temperatures, the 
overall hopping potential contribution to the total action is relevant also at low T. Precisely this property causes the 
anomalous upturn in the heat capacity linear coefficient. Summing over a large number of paths is essential to recover 
the correct thermodynamical behavior of the heat capacity in the zero temperature limit. 

Further investigation also reveals that the upturn persists both in the extremely narrow (J ~ lOmeV) and in the 
wide band (J ~ leV) regimes. Our method accounts for a variable range hopping on the r scale which corresponds 
physically to introduce some degree of disorder along the linear chain. This feature makes our model more general 
than the standard SSH Hamiltonian with only real space nearest neighbor hops. Although I am not aware of any other 
direct computation of specific heat in the SSH model, hopping type mechanisms have been suggested [31] to explain 
the striking conducting properties of doped polyacetylene at low temperatures. Since the specific heat directly probes 
the density of states and integrating over T the specific heat over T ratio one can have access to the experimental 
entropy, the method here presented may provide a new approach to analyse the transition to a disordered state which 
indeed exists in polymers [32]. In this regard it is worth remarking that glassy systems [33] infact exhibit a low T 
upturn in the specific heat over T ratio due to tunneling states for groups of atoms providing a non magnetic internal 
degree of freedom in the potential structure [34,35] . 

5. CONCLUSIONS 

We have developed a general path integral method to study the interplay between electron motion and oscillators 
bath in a semiclassical Su-Schrieffer-Hecger model. Using the fact that the time dependent average energy is linear in 
the oscillator displacement, we have analytically deduced the full partition function as a functional of the e-ph source 
term and computed the two particle correlations induced by the phonon baths at some selected temperatures: at 
increasing temperatures the electrons are correlated over larger distances. Then, we have focusscd on the perturbing 
effect of the electron particle motion on the thermodynamics of a single oscillator. After expanding the time dependent 
oscillator path in Fourier series, we have derived the k— order cumulant correction (eq.(lO)) to the oscillator partition 
function in the case of a time averaged electron particle path. At decreasing temperatures an increasing number (Np) 
of Fourier components in the oscillator path is required to compute the cumulant terms with accuracy. Np is also a 
growing function of the oscillator energy. The proposed method permits to evaluate, for a given e-ph coupling, the 
cumulant corrections to the harmonic partition function and free energy of each oscillator at any temperature. While 
higher k— order terms become more relevant at low T and they have to be included in the partition function series 
expansion, the main corrections to the harmonic free energy arc instead ascribable to the lower fc-ordcr cumulants and 
their effect is also evident at increasing temperatures. Further refinements on the presented results may be obtained 
by integrating the oscillator partition function over an ensamble of r dependent particle paths. Finally, the path 
integral method has been applied to compute the heat capacity of the system as a function of the e-ph coupling and 
the oscillators energies. We find a peculiar upturn in the low temperature plots of the heat capacity over temperature 
ratio indicating that a glassy like behavior can arise in the linear chain as a consequence of a time dependent electronic 
hopping with variable range. 
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FIG. 1. Square root of the two-points correlation functions in units A. A bath of ten phonon oscillators is considered, the 
largest phonon energy being hujio — AOmeV (intermediate phonon spectrum) and fiuJio = GOmeV (high phonon spectrum), 
(a) T= IK, (b) T= lOOK, (c) T= 300K. 



FIG. 2. Square root of the two-points correlation functions in units A plotted, at different temperatures, in the range 
T < 0.04meV~^ . An intermediate phonon spectrum is considered, the largest oscillator energy being hojio = 4QmeV . 



FIG. 3. (a) Neperian logarithms of the first three cumulants in oq.(fO). (b) Phonon free energy and anharmonic corrections 
due to the first three cumulants. The e-ph coupling a is in units eVA~^. w is the phonon energy. 



FIG. 4. As in Fig.3 but with a large e-ph coupling a. 
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FIG. 5. Total heat capacity over temperature for three vahies of the effective couphng x (hi units meV^). A bath of ten 
phonon oscillators is considered, the largest phonon energy being huiio = 20meV. The phonon heat capacity is also plotted. 



FIG. 6. Total heat capacity over temperature for five values of the eflfective coupling x (in units meV^). The largest oscillator 
energy of the phonon bath is huiio = 4QmeV . 
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